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Abstract 

We investigate nuclear matter properties in the relativistic Brueckner approach. The in- 
medium on-shell T-matrix is represented covariantly by five Lorentz invariant amplitudes from 
which we deduce directly the nucleon self-energy. To enforce correct Hartree-Fock results we 
develop a subtraction scheme which treats the bare nucleon-nucleon potential exactly in accor- 
dance to the different types of meson exchanges. For the higher order correlations we employ two 
different covariant representations in order to study the uncertainty inherent in the approach. 
The nuclear matter bulk properties are only slightly sensitive on the explicit representation used. 
However, we obtain new Coester lines for the various Bonn potentials which are shifted towards 
the empirical region of saturation. 

Introduction 



he investigation of nuclear matter properties within the relativistic Dirac-Brueckner-Hartree-Fock 
(DBHF) approach [I], 0, |3|, [|, ^ |5| remains a fundamental topic in theoretical nuclear structure studies. 
Compared to non-relativistic approaches the relativistic DBHF treatment turned out to be a major 
step forward in the explanation of the saturation mechanism of nuclear matter. The saturation points 
obtained for non-relativistic calculations, throughout all possible choices of different nucleon-nucleon 
interactions, are located on the so called 'Coester line' @ which does not meet the empirical saturation 
region. Using modern nucleon-nucleon interactions of the one-boson exchange type |8j the relativistic 
calculations also reveal such Coester lines which are, however, significantly shifted towards the empirical 
region ||. 

On the other hand, many details of the relativistic theory are still not fully resolved. In particular, 

the precise form of the nucleon self-energy, i.e. the magnitude and the momentum dependence of the 

scalar and vector self-energy components are a question of current debate ||, |], |5| . Since the self-energy 

describes the dressing of the particles inside the medium and thus determines the relativistic mean 
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field this fact states a severe problem. Different techniques to handle the DBHF problem can lead to 
significantly different results 0, [|, |5|, ^| . In a recent work H we found that the momentum dependence 
of the nucleon self-energy is dominated by the one-pion exchange contribution which accounts for the 
nuclear tensor force. 

Unfortunately, the treatment of the 7rNN vertex and the corresponding self-energy contributions is 
closely connected to a severe ambiguity in the T- matrix representation || . The DBHF approach starts 
from a realistic nucleon-nucleon potential of the one-boson exchange type, i.e. the Bonn potentials 
0. As for the free two-body scattering problem, anti-particle states are neglected, and thus one 
works exclusively with positive energy states. Hence, a direct determination of the nucleon self- 
energy operator is not possible since not all matrix elements of this operator are known. Horowitz 
and Serot have therefore developed a projection technique to determine the scalar and vector self- 
energy components from the in-medium T- matrix [|TJ. In this approach the T- matrix is represented 
covariantly by Dirac operators and Lorentz invariant amplitudes where the latter are determined from 
the positive-energy on-shell T-matrix elements. 

The whole problem arises from the fact mentioned above, namely that one does not include 
negative energy states and therefore neglects the excitation of anti-nucleons. The inclusion of negative 
energy excitations with 4 states for each spinor yields 4 4 = 256 types of two-body matrix elements 
concerning their spinor structure. Symmetry arguments reduce this to 44 for on-shell particles. ]T0 |. 
If one takes now only positive energy solutions into account this reduces to 2 4 = 16 two-body matrix 
elements. Considering in addition only on-shell matrix elements the number of independent matrix 
elements can be further reduced by symmetry arguments down to 5. Thus, all on-shell two-body 
matrix elements can be expanded into five Lorentz invariants. But these five invariants are not unique 
since the Dirac matrices involve always also negative energy states and thus a decomposition of the 
one-body nucleon-nucleon potential into a Lorentz scalar and a Lorentz vector contribution depends 
on the choice of these five Lorentz invariants mentioned above. The best choice would be to separate 
completely the negative energy Dirac states. But since this is not possible, there is not a unique 
but only an 'optimal choice'. The topic of this paper is the form of this 'optimal choice' of the five 
invariants. 

The paper is organized as follows: In section 2 we briefly review the Dirac-Brueckner Hartree- 
Fock approach. In section 3 we introduce the projection technique and discuss the different covariant 
representations used for the on-shell T-matrix. Nuclear matter bulk properties are then discussed in 
section 4. At the end we summarize and conclude our work. 



2 The relativistic Brueckner approach 
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2.1 The coupled set of equations 

In the relativistic Brueckner approach the nucleoli inside the nuclear medium is viewed as a dressed 
particle in consequence of its two-body interaction with the surrounding nucleons. The in-medium 
interaction of the nucleons is treated in the ladder approximation of the relativistic Bethe-Salpeter 
equation 

T = V + i J VQGGT , (1) 

where T denotes the T-matrix. V is the bare nucleon-nucleon interaction. The intermediate off- 
shell nucleons in the scattering equation are described by a two-particle propagator iGG. We replace 



this propagator by the Thompson propagator [TTJ. The Pauli operator Q in the Thompson equation 



accounts for the influence of the medium by the Pauli-principle and projects the intermediate scattering 
states out of the Fermi sea. 

The Green's function G fulfills the Dyson equation 

G = G + G m . (2) 

Gq denotes the free nucleon propagator while the influence of the surrounding nucleons is expressed 
by the nucleon self-energy S. In Brueckner theory this self-energy is determined by summing up the 
interaction with all the nucleons inside the Fermi sea 

S = -i J (Tr[GT] - GT) . (3) 

F 

The Dirac structure of the self-energy in isospin saturated nuclear matter follows from translational 



and rotational invariance, parity conservation and time reversal invariance fll^] . In the nuclear matter 
rest frame the self-energy has the simple form 

E(k, k F ) = E 8 (k, k F ) -7 S o (A;, k F ) +7-kS v (A;, k F ) , (4) 

with kfj, being the nucleon four-momentum. By taking the traces in Dirac space as |J] 

X B = \tr\E] , £ o = ^r[ 70 X] , S v = ^tr[ 7 -kS] (5) 

one can calculate the different Lorentz components of the self-energy. 



2.2 The in-medium T-matrix 



We apply the relativistic Thompson equation |LT| to solve the scattering problem of two nucleons in 



the nuclear medium. In the two-particle center of mass (cm.) frame - the natural frame for studying 
the two-particle scattering process - this Thompson equation can be written as @. |J 

T(p,q,x)\ c . m . = V(p,q) (6) 

J (2vr) 3 ' E* 2 (k)2E*{q)-2E*{k)+te V ; ' 
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where q = (qi — c^)/2 is the relative three-momentum of the initial state while k and p are the 
relative three-momenta of the intermediate and final states, respectively. The total four-momentum of 
the two-nucleon system is P* = q\ + q%. Vs* — 2E*(q) = 2^q 2 + m*p is the starting energy in (^]). If 
qx and q2 are nuclear matter rest frame momenta of the nucleons in the initial state, the boost- velocity 
u into the cm. frame is given by 

u = P/v^ + P 2 , (7) 

with the total three- momentum and the invariant mass P = qi + q2 and s* = (E*(qi) +E*(q^2)) 2 — P 2 , 
respectively. In Eq. ([?]) x denotes the set of additional parameters x = {kp,rh* F , |u|} on which the 
T-matrix depends. 

Applying standard techniques as explained in detail by Erkelenz || we solve the Thompson 
equation in the cm. frame and calculate the plane-wave helicity matrix elements of the T-matrix. 
The subspace of negative energy states is omitted in the current Brueckner approach. In this way we 
avoid the delicate problem of infinities in the theory which would generally appear if we would include 
contributions from negative energy nucleons in the Dirac sea jl], [(J . 



3 Covariant representations and the nucleon self-energy 
3.1 Pseudo-scalar representation 

To use the trace formulas, Eqs. (^), one has to represent the T-matrix covariantly. A set of five 
linearly independent covariants is sufficient for such a T-matrix representation because on-shell only 
five helicity matrix elements appear as solution of the Thompson equation. A linearly independent 
although not unique set of five covariants is given by the Fermi covariants 

S = 1 ® 1, V = 7^ ® 7m , T = ® <v, A = 757^ ® 757^, P = 75 ® 75- (8) 

Using this special set - the so called 'pseudo-scalar choice' - the on-shell T-matrix for definite isospin 
I can be represented covariantly as [[jj 

T\\p\,9,x) = F s I (|p|,^,x)S + F|(|p|,^,x)V + Fi(|p|,^,x)T 

+ Fi(|p|,0,oOA + F£(|p|,M)P • (9) 

Here p and 9 denote the relative three-momentum and the scattering angle between the scattered 
nucleons in the cm. frame, respectively. Applying the covariant representation (|S|) for the on-shell 
T-matrix the nucleon self-energy in isospin saturated nuclear matter is evaluated to be 0] 

v a i ^ f d3( * e ( kp ~ 1*11) f~* i Ei, x* z? 1 nn\ 
Z a p{k,k F ) = J ( 27r ) 3 E*(q) [ m *' la f }Fa+ fi*P F v\ ' ( 10 ) 
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where the isospin averaged amplitudes are defined as 



Fi(|p|,0,z) := - [^ I=0 (|p|,0,x) + 3^ I=1 (|p|,0,x) 



(11) 



In Fig. [l| we show the result of a self-consistent DBHF calculation for the nucleon self-energy in 
nuclear matter applying as representation for the on-shell T- matrix the ps representation (|9[). As 
bare interaction we have used the Bonn A potential || and, for comparison, the a-u model potential 
which was originally used by Horowitz and Serot [0J. The Fermi momentum is kp = 1.34/m" 1 . As 
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Figure 1: Momentum dependence of the DBHF nucleon self-energy components in nuclear matter at 
kp = 1.34/m -1 using as bare nucleon-nucleon potential Bonn A (solid) and the a-u model potential 
(dashed). For the T-matrix the ps representation (|9|) is applied. 

already discussed in Ref. Q, we see a pronounced momentum dependence of the nucleon self-energy 
components with the full Bonn A while in the case of the a-uo model potential the dependence on the 
momentum is rather weak. A strong momentum dependence leads to unphysical results deep inside 
the Fermi sea since the effective mass drops to values which are close to zero. Therefore in Ref. ]5| 
the strong momentum dependence of the self-energy was studied in more detail and found to originate 
mainly from the one-pion exchange contribution to the self-energy. 



3.2 Complete pseudo-vector representation 

To suppress the undesirable pseudo-scalar contribution of the pion to the nucleon self-energy we have 
to use the 'complete' pv representation of the T-matrix [13 

T^lpl^x) = ^(|p|,^,x)S-4(|p|,^,x)S + ^(|p|,^,x)(A-A) 

+ ^v(|p|,0,a;)PV-^(|p|,0,a;)PV . (12) 

The amplitudes g l (0) are explicitly given in |14]]. In Fig. |] we present the full self-consistent DBHF 



calculation with the 'compete' pv representation of the T-matrix ||. The DBHF nucleon self-energy 
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Figure 2: Momentum dependence of the DBHF nucleon self-energy components in nuclear mat- 
ter at kp = 1.34/m" 1 using Bonn A as bare nucleon-nucleon interaction. For the T-matrix the ps 
representation (Eq. (Bf), solid) and the 'complete' pv representation (Eq. (|T2"D, dashed) are applied. 

components are indeed weakly momentum dependent. The single-pion exchange contribution to the 
interaction, which was previously dominating at low nucleon momenta, is now strongly suppressed. 
Consequently the result within the 'complete' pv representation using the Bonn A potential resembles 
the result within the a-uj model potential, see Fig. [I] where the ps representation was used. To suppress 
the pion contribution to the in-medium T-matrix a correct pseudo-vector like covariant representation 
is essential for the calculation of the nucleon self-energy in nuclear matter. 




3.3 Covariant representations of the subtracted T-matrix 

The 'complete' pv representation successfully reproduces the HF nucleon self-energy in the case of 
the pion exchange. However, as already pointed out in ||, the 'complete' pv representation fails to 
reproduce the HF nucleon self-energy if other meson exchange potentials are applied as bare interaction. 
Hence, it should be reasonable to treat the bare interaction and the higher order ladder graphs of the 
meson exchange potential separately. Since the single-meson exchange potential is actually known 
analytically we can represent it covariantly by a mixed representation of the form 

V = V^ + V a p ^ s . (13) 

Here the 7T- and //-meson contributions are treated as pseudo-vector while for the (a, u, p, 5)-meson 
contributions of the Bonn potential the ps representation is applied. The higher order ladder diagrams 
of the T-matrix 

Tsub = T — V = i J VQGGT = £ / V(iQGGV) n , (14) 

n=l 
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in the following called the subtracted T-matrix, can not be represented correctly in a mixed form 
since we can not disentangle the different meson contributions to this part of the full in-medium 
interaction. The representation of the subtracted T-matrix remains therefore ambiguous. However, 
if the pion exchange dominantly contributes to the Hartree-Fock level a ps representation of the 
subtracted T-matrix should be more appropriate because then the higher order contributions of other 
meson exchange potentials are not treated incorrectly as pseudo-vector. Thus the most favorable 
representation of the T-matrix is given by the ps representation 

T P = T! ub + V^ + V a p ^ s . (15) 

Here the ps representation for Tg ub is determined via the matrix elements 

< P A' 1 A;|Tl ub (x)|qA 1 A 2 >:=< pX'X\T J {x) - ^(aOlqAiAa > , (16) 

with subsequently applying the projection scheme as in Eq. (??). An alternative representation of the 
T-matrix is given by the pv representation 

T PY = TZ + V^ + V^ S , (17) 

where the subtracted T-matrix is represented by the 'complete' pv representation fll2l) . This represen- 
tation is similar to the 'complete' pv representation of the full T-matrix, however, with the advantage 
that now the pseudo-scalar contributions in the bare nucleon-nucleon interaction, e.g. the single-omega 
exchange potential, are represented correctly. In the next section we will use both representations, 
dl5|) and (17), to study the properties of nuclear matter in the DBHF approach. In this way we can 



determine the influence of the higher order ladder graphs to the in-medium interaction in a more 
quantitative way. Furthermore, these two representations set the range of the remaining ambiguity 
concerning the representation of the T-matrix, i.e. after separating the leading order contributions. 



4 The equation-of-state of nuclear matter 

In the relativistic Brueckner theory the energy per particle is defined as the kinetic plus half the 
potential energy 

1 1 m*(k) 

£M = -E<«A(k)| 7 -k + M+-£(AOK(k)>^-M . (18) 

In Fig. |3] we show the binding energy per particle E/A as a function of the density, calculated with 

Bonn A, B and C. For the T-matrix the subtraction scheme with the ps representation ( P~5|) is applied. 

With Bonn A one can reproduce the empirical saturation point of nuclear matter, shown as shaded 

region in the figure. The other Bonn potentials give less binding energy although the saturation density 
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Figure 3: Binding energy per particle as a function of nuclear matter density. As bare nucleon-nucleon 
interaction the potentials Bonn A, B and C are used. For the T-matrix the subtraction scheme with the 
ps representation QT5| ) is applied. The shaded box denotes the empirical saturation region of nuclear 
matter. 

is always close to the empirically known value. The result for the binding energy per particle using the 
two representations (|15|) and flT7|) for the T-matrix are very similar. At saturation density the binding 
energy is only 0.5 MeV smaller using the pseudo-vector representation of the subtracted T-matrix. 
Thus, the energy per particle is not very sensitive on the explicit representation of the subtracted 
T-matrix. 

The present results are summarized in Fig. [| where the corresponding saturation points for the 
three different versions of the Bonn potential are shown. We compare the results with the two repre- 
sentation of the subtracted T-matrix with the results of the calculation of Brockmann and Machleidt 
(BM), Ref. p|. With the improved representation schemes ([15]) and ( |TTD for the T-matrix one obtains 
new 'Coester-lines' which are left of the original one, i.e. shifted towards the empirical region. The 
refined treatment of the T-matrix representation leads to an enhancement of the binding energy con- 
nected with a reduced saturation density. As in the previous calculations, Bonn A is still the only one 
which meets the empirical region. However, due to an increased binding energy Bonn B is now much 
closer to empirical region. This observation is consistent with the present treatment. The different 
types of Bonn potentials essentially vary in the strength of the nuclear tensor force determined by the 
ttNN form factor. Bonn A which has the smallest tensor force yields the smallest D-state probability 
of the deuteron and only a pure description of the 3 Di phase shift 0, || . Thus it appears that a refined 
treatment of the pion exchange leads to improved nuclear matter results for the more realistic Bonn 
B potential. Bonn C, however, is still far off the empirical region. 

Furthermore, it can be seen from Fig. (^) that the final nuclear matter bulk properties depend 

8 



-10 



1.2 1.3 1.4 1.5 

k F [fm 1 ] 

Figure 4: Saturation points of nuclear matter. As bare nucleon-nucleon interaction the Bonn poten- 
tials A,B and C are used. For the T-matrix the subtraction scheme with the ps representation (Eq. 
(|ToD, filled circles) and the pv representation (Eq. (|TTD, open circles) are applied. As open triangles 
the results of the calculation of Brockmann and Machleidt (BM) , Ref . || , are shown. 

only moderate on the representation of the subtracted T-matrix. In Ref. || we tried already to 
determine the range of inherent uncertainty in the relativistic Brueckner approach which is due to 
the ambiguities concerning the representation of the T-matrix discussed in Section 3. By the separate 
treatment of the Born contribution to the T-matrix we end up now with a much narrower uncertainty 
band which is given by the ps or complete pv representation of the ladder kernel, i.e. the subtracted 
T-matrix. Over the different types of Bonn interactions the two methods lead to a variation of 0.5 
MeV in the binding energy, 0.1-0.2 fm" 1 in the Fermi momentum, and to about 30 MeV concerning 
the value of the effective mass at saturation density. The values for incompressibility are also close 
in the two approaches, i.e. they differ by less than 10 MeV. Within the ps representation of the 
subtracted T-matrix, Bonn B and C now yield very small kompression moduli around K = 150MeV 
and K = 115MeV, respectively. For Bonn A a value of K = 230MeV is obtained. This value agrees 
with the empirical value of the kompression modulus of K = 210 ± 30MeV. Here Brockmann and 
Machleidt found much larger values for all three Bonn potentials. 
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5 Summary 

We have investigated the nuclear matter properties in the relativistic Brueckner approach. The re- 
quired representation of the T-matrix by Lorentz invariant amplitudes suffers thereby from on-shell 
ambiguities concerning the pseudo-scalar or pseudo-vector nature of the interaction. We minimized 

this ambiguity by separating the leading order, i.e. the single-meson exchange, from the full T-matrix. 
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The remaining higher order correlations, i.e. the ladder kernel, are then represented either completely 
as pseudo-scalar or as pseudo- vector. 

As a major result of our investigation we obtain new 'Coester lines' for the various Bonn po- 
tentials. Compared to previous treatments these are shifted towards the empirically know saturation 
point. Bonn A is still the only potential which really meets the empirical region of saturation, but, 
with improved saturation properties compared to previous treatments. The refined treatment of the 
pion exchange leads on the other hand also to improved results for the - from the view of the phase 
shift analysis - more realistic Bonn B potential. Furthermore, we found that the equation-of-state is 
strongly softened compared to previous calculations. Actually with Bonn A we obtain a kompression 
modulus of K ~ 230MeV which is in good agreement with the empirical value. 

To summarize our results, we obtained new results for the nuclear matter properties within the 
projection technique employing an new method for the T-matrix representation. The final results are 
at lower densities almost insensitive on the explicit choice made for the representation. However, at 
higher densities, certain differences occur when using different representation schemes. We want to 
stress that the ambiguity in the projection technique is still not fully resolved yet. We plan to look 
on off-shell T-matrix elements in the future since off-shell matrix elements of the pseudo-scalar and 
pseudo-vector covariants differ significantly. We hope that this might bring more insight on what is 
the correct on-shell representation of the T-matrix. 
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